Does M31 result from an ancient major merger? 
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ABSTRACT 

The M31 haunted halo is likely associated with a rich merger history, currently assumed to be 
caused by multiple minor mergers. Here we use the GADGET2 simulation code to test whether 
M31 could have experienced a major merger in its past history. Our results indicate that a 
(3±0.5):1 gaseous rich merger with r per i ceri ter=25±5 kpc and a polar orbit can explain many 
properties of M31 and of its halo. The interaction and fusion may have begun 8.75±0.35 and 
5.5 ±0.5 Gyr ago, respectively. Observed fractions of bulge, and the thin and thick disks can be 
retrieved for a star formation history that is almost quiescent before the fusion. This also accords 
well the observed relative fractions of intermediate age and old stars in both the thick disk and 
the Giant Stream. In this model, the Giant Stream is caused by returning stars from a tidal tail 
which contains material previously stripped from the satellite prior to the fusion. These returning 
stars are trapped into elliptical orbits or loops for long periods of time which can reach a Hubble 
time, and are belonging to a plane that is 45 degrees offset from the M31 disk position angle. 
Because these streams of stars are permanently fed by new infalling stars with high energy from 
the tidal tail, we predict large loops which scale rather well with the features recently discovered 
in the M31 outskirts. 

We demonstrate that a single merger could explain first-order (intensity and size), morpho- 
logical and kinematical properties of the disk, thick disk, bulge and streams in the halo of M31, 
as well as the observed distribution of stellar ages, and perhaps metallicities. This challenges the 
current scenarios assuming that each feature in the disk (the 10 kpc ring) or in its outskirts (thick 
disk, the Giant Stream, and the numerous streams) is associated with an equivalent number of 
minor mergers. Given the large number of parameters, further constraints are certainly required 
to better render the complexity of M31 and of the substructures within its halo which may ulti- 
mately lead to a more precise geometry of the encounter. This would allow us, in principle, to 
evaluate the impact of such a major event on the Andromeda system and Local Group. 
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1. Introduction 

Our nearest giant neighbour, M31 has attracted 
a considerable amount of interest since the discov- 
ery of many large-scale structures that surround 



its outskirts. These prominent structures include 
the Giant Stream (infall of 1.5 x 10 s M Q stel- 



lar mass, Ibata et al.||2001 ) and the gigantic thick 
disk containing about 10% of the disk luminosity 
(llbata et al.||2005b. The halo of M31 is haunted 
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by not less than 16 substructures that led |Tanaka| 



et al. (20101 to conjecture that they are due to as 



many accretion events involving dwarf satellites 
with mass 10 7 -10 9 M ffl since z ~ 1. Observations 



at radio wavelengths (Westmeier et al. (20051; 
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Thilker et al. (2004); and references therein) reveal 



high velocity clouds, whose location and kinemat- 
ics partly follow the Giant Stream. Moreover, a 
lOkpc, pseudo-ring inserted in the M31 disk dom- 
inates star formation (Baade & Arp 1964), HI 

gas 



2006 



Roberts 1966), molecular gas (Nieten et al. 
and dust emission dGordon et al.||2006[ ). 



Simulations of the numerous structures in the 
M31 outskirts have always assumed them to be 
caused by minor satellites. M32 is understood to 



be the perturber of the spiral arms (Byrd 1978 



1983), while NGC 205 has been modelled as the 



cause of the warp seen in the optical and HI disks 
of M31 dSato & Sawa 1119861). The thick disk has 



been modelled by earlier disruption, several Gyr 
ago, of dwarf galaxies on prograde orbits that are 
coplanar with the disk (iPenarrubia et al. 2006). 



Simulations of the Giant Stream assumed a col- 
lision with an unknown satellite ~ 0.7 Gyr ago 
( [Font et al.||2006| |Fardal et al"]|2008| . The 10 kpc 
ring could also be reproduced by an interaction 



with M32, assuming a polar orbit (Block et al. 



2006) 



In spite of their success in reproducing M31 
structures, several of these simulations may be 
speculative because the M31 satellite orbits are 
currently unknown (e.g., Fardal et al.|2009 |. More 
problematic is the fact that stars in the Giant 



Stream (Brown et al. 2007) have ages from 5.5 



to 13 Gyr, which is difficult to reconcile with a 
recent collision (e.g., Font et al.|2008 ). If the sub- 
structures are formed from different progenitors, 
why do they show obvious similarities in metallic- 
ity (e.g., Ferguson et al. 2005) ? The main moti- 



vation for only investigating minor mergers is to 
preserve the M31 disk (iMori & Rich||2008|) . But 



the M31 disk is not necessarily very old and per- 
manent, as its associated stellar clusters display 



young to moderate ages (< 5-7 Gyr, Beasley et 



al. 2004), and most of the stars in the outskirt 
structures are older than this. One should not 
discard the possibility that the gigantic structures 
may have been formed at the same time or even 
earlier than the disk. In fact many co-workers in 



the field (Rich 2004 Ibata et al. 12004 Brown et al 



2006 ) have hypothesised a possible major merger 
in the past history of M3]Q More recently, the 



1 Quoting van den Bergh ] < |2005| >: ' 'Both the high metal 



licity of the M31 halo, and the r 1 / 4 luminosity profile of 



kinematics of the M31 globular system has been 
attributed to an ancient major merger (e.g 

polo] ). 



Bekki 



Understanding the nature of M31 certainly has 
an impact in cosmology because it and the Milky 
Way, are the only two giant spirals in our immedi- 
ate neighbourhood. In contrast to the Milky Way, 
M31 has properties (absolute K luminosity, angu- 
lar momentum, [Fe/H] in the outskirts) similar to 
the average of large spirals having the same ro- 
tation velocity ( Hammer et al.p007 ). Half of the 
progenitors of spiral galaxies in this range were not 
relaxed 6 billion years ago, according to the de- 
tailed studies of their morphologies (IHammer et 



al. 



al. 



[2005] |Neichel et aTJ|2008{ |Delgado-Serrano~et 



2010) and kinematics (Puech et al.[2008 Yang 



et al. 2008). Hammer et al. (2009) verified that 



various merger phases can reproduce quite well 
the observed morphologies and kinematics of these 
non-relaxed galaxies. This agrees with the disk re- 
building scenario ( |Hammer et al.|2005[ ), according 
to which most spiral disks have been rebuilt follow- 
ing a major, gas-rich merger during the past 8-9 
Gyr. High gas fractions have been shown to be 



essential during this process (Hopkins et al. 2008 



2010). Understanding the different substructures 



of M31 as resulting from a single event such as a 
gas-rich major merger is an important step in val- 
idating or dismissing such a channel for spiral disk 
formation. 

The goal of this study is to investigate whether 
a past major merger, instead of multiple minor 
mergers, can reproduce most of the peculiarities 
of M31. The ACDM cosmology {H =70 km s" 1 
Mpc" 1 , flM = 0.3 and = 0.7) is adopted 
throughout the paper. 

2. Observational constraints and initial 
conditions 

We assume in the following that M31 results 
from a past major merger. Our goal is to simul- 
taneously reproduce: (1) the M31 disk at its rota- 
tion velocity, and the bulge, with B/T=0.28; (2) 
the presence of a very extended, co-rotating thick 
disk; (3) the lOkpc dust, HI, star-forming ring; (4) 



the Andromeda galaxy, suggest that this object might have 
formed from the early merger and subsequent violent relax- 
ation, of two (or more) relatively massive metal-rich ances- 
tral objects." 
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Fig. 1. — Chronological sketch of the structures surrounding M31. In the central panel (reproduced from 
Ibata et al.|2005 1, the large and thick rotating disk is a vast flattened structure with a major axis of about 



4 degrees. Squares represent fields observed by Brown et al. (|2006 
measurements by arrows 



2007 2008), and are linked to their 



the Giant Stream, and (5) the observed age dis- 
tributions of stars in the different substructures. 
Such a choice may appear somewhat subjective 
although it has been adopted after thoroughly re- 
viewing many of the M31 properties. First, be- 
cause the M31 system is very complex to model, 
we need to limit the number of features that can be 
reproduced. Second, the limitations of our model 
(number of particles) prevent us from describing 
very small structures, for example the double nu- 
cleus of M31 (but see Hopkins fc Quataert||2010 |. 
Third we aim at reproducing the zero and first- 
order features of each substructure of M31, includ- 
ing their stellar age, mass, and kinematics (if they 
are available from observations), but not the de- 
tailed morphology of each of them (see discussion, 
section 5). 



According to Hopkins et al. (2010) (see their 
Figure [7]), a mass ratio from 0.3 to 0.5 is required 
to reform a new bulge with B/T^0.3 in the rem- 



nant galaxy for gas fraction^] ranging between 
0% and 50%. However, reforming a significant 
thin disk requires that enough gas be preserved 
after the fusion of the cores: it is mostly made 
of stars from the gaseous thin disk that immedi- 



ately reforms after the collision (Barnes, J.E. 2002 



Abadi et al. 


2003 


kins et al. 


2009 


1. 



Governato et al.||2007| |Hop-| 



Thus reforming an Sb galaxy 
like M31 probably requires a gas fraction in ex- 
cess of 50% in the progenitors. The total bary- 
onic mass of M31 is l.lxlO 11 M (Hammer et al. 



2007), and at first approximation we can assume 
that this is the value for the sum of the two pro- 
genitor masses. In order to prevent a too violent 
relaxation at the center of the main progenitor, 
we adopt a prograde-retrograde orientation of the 
spin axes, which has been shown to be more favor- 



2 The gas fract ion in the progenitors should reach a higher 
value than in |Hopkins et al.| l |2010| | , since they define the 
gas fraction as its value just before the fusion. 
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able in rebuilding a thin disk after a major merger 
( |Hopkins et~aL|2008| ). 

The presence of a gaseous ring favours a polar 
orbit. As shown bylMori & Richl (|2008|), the Giant 



Stream may result from the returning material of 
the tidal tail formed just before the last passage 
of the satellite. A similar although longer-lived 
tidal tail is predicted for a 3:1 major merger with 
a polar orbit, which is associated with the pas- 
sage of the secondary galaxy just before fusion. 
During the first passage and until the fusion, star 
formation is especially enhanced in the secondary 



galaxy that is harassed by the main one (e.g., Cox 



et al. 2008). If the Giant Stream was associated 



with returning particles from the tidal tail formed 
during the second passage (near fusion), it would 
contain mostly stars older than the epoch of the 
fusion: this is because the star formation cannot 
hold for a long time within tidal tails due to their 
expected dilution. 

We thus propose from Fig. [I] a chronological 
history of the different structures in M31, as this 
figure can be used as a clock for determining the 
occurrence of merger phases. The star formation 
history of the whole system is enhanced during the 
first passage until the fusion and then at the fu- 
sion itself (see Cox et al. 2008). During a merger 
event, most of the gas and stars in the remnant 
outskirts have been deposited by tidal tails formed 
during the first passage and later, during the fu- 
sion of the cores. A few hundred billion years af- 
ter its formation, the tidal tail dilutes, provoking 
a natural quenching of the residual star formation 
(see Wetzstein et al. 2007). Thus the age of the 
material brought by tidal tails provides, with a 
relatively small delay, the date of both the first 
passage and fusion times. In the following, we as- 
sume that the first passage occurred from 8.5 to 
9 Gyr ago, and that the corresponding tidal tails 
are responsible for the halo enrichment seen in the 
21 kpc and 35 kpc fields, without significant star 
formation more recent than 8.5 Gyr. The thick 
disk has a star formation history comparable to 
that of the Giant Stream, and is also generated 
by material returning to the galaxy mostly from 
tidal tails generated at the fusion. Because their 
youngest significant population of stars has ages of 
5.5 Gyr, the delay between the first passage and 
fusion ranges between 3 and 3.5 Gyr. This could 
be accommodated for by relatively large impact 



parameters (20 to 30 kpc). 

From the above, we can settle the initial con- 
ditions for a major merger assumed to be at the 
origin of M31. Table[T]describes the adopted phys- 
ical parameters of such an interaction. Given the 
huge amount of observational data, it is beyond 
the scope of this study to reproduce the details 
of all the numerous structures in the M31 system. 
Instead, our aim in using hydrodynamic simula- 
tions is to determine whether or not these numer- 
ous substructures can be attributed to a single ma- 
jor merger in the past history of M31. 

3. Simulations 

We used the GADGET2 hydrodynamical code 



(Springel 2005) and initial conditions similar to 
that of |Cox et all (120061 120081) . For the dark mat- 



ter, we adopted a core density profile as in |Barnes^1 
(2002), with a core size of 5.3 kpc and 3.06 



J.E 



kpc for the main galaxy and the satellite, respec- 
tively. We verified that our results are not signif- 
icantly affected by changing the density profile to 
a Hernquist model that is quite similar in central 
regions to a Navarro-Frenk- White model (see sec- 
tion 5). Concerning the implementation of feed- 
back and cooling, we have verified step by step our 
ability to reproduce both the isolated and merger 
cases for all the different combinations of feedback 



and cooling in Cox et al. (2008), and as such, all 



parameters used in our simulations are very sim- 
ilar to those of Cox et al. (see Wang 2010 in 



preparation) Usually, the free parameters describ- 
ing the star formation efficiency and the feedback 
strength are chosen in order to match the Schmidt- 
Kennicutt law between star formation and gas sur- 
face densities (see Cox et al.|2006 ). Given the rel- 
atively large scatter of this relation, a large num- 



ber of combinations can be accommodated (Cox 



et al. 2006 2008). In the following, we explored 



several combinations that are in agreement with 
this relation, but with the additional constrain of 
preserving a significant gas reservoir in the pro- 
genitors, before fusion (see Sect. 4.2). Properties 
of the progenitors are listed in Table [l] and these 
galaxies have been generated to follow the bary- 
onic Tully Fisher relation (see Puech et al. 2010 



for an argument in favour of a non-evolving rela- 
tion) . Their atomic gas content has been assumed 
to be three times more extended than the stars, as 
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Fig. 2. — Different phases of a 3:1 major merger for M31 (r per i center =24 kpc, Gall incy=70, Gal2 incy=-110). 
The simulation starts 9.3 Gyr ago (T=0, z=1.5) and the first passage occurs 0.7 Gyr later. Then new stars 
are forming (green color), especially in the secondary galaxy until the fusion, which occurs at 4.5 Gyr. The 
elapsed time between the first passage and fusion is 3.8 Gyr, as the pericenter radius is large. During the 
second passage (T=4.2 Gyr) a tidal tail containing many newly formed, intermediate-age stars (green dots) 
is formed. Later on this material returns to the galaxy forming the Giant Stream, enriched with stars formed 
from 5 to 8 Gyr ago. The resulting galaxy in the last panel (T=9.3 Gyr) is compared with the inserted 
M31 image at the same scale (from Ibata et al.||2005 , see an enlarged view of this insert on bottom- right). 
The Figure also illustrates the formation of the Giant Stream (see section 4.4) in a case for which the two 
galaxies have angle difference near the resonance (180 degrees) providing many particles stripped from the 
satellite at 4.2 Gyr. Some of the particles within the tidal tail have velocity below the escape velocity of 
the remnant system and are gradually falling back to the galaxy. They are tracing loops around the newly 
formed disk as indicated by the arrow in the panel at T=6.8 Gyr. Loops are fed by newly-coming particles 
falling back from the tidal tail and are persistent until 9.3 Gyr and later. The dotted rectangle in the 4.2 
Gyr panel illustrates how we have selected the tidal tail particles (see section 4.4). 
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Table 1: Initial and Adopted Conditions for a Major Merger Model for M31. 



Ingredient 


xveoieQ xvian^e 


Comments 


rvQUpLCQ Ivail^e 


total mass 


^ ^ v 1 nil A/f 

0.0 X 1U -i*^0 


20% of baryons 




vn a o o voii f~\ 
llldoo IdLlU 


9 A 


tn rpfnrm R /T^D "? 


9 ^ 3 ^ 


f Call 
J gas vjcUI 




p"VT"ipptprl fit 7 — 1 ^ a 

CAJJC^lCd at Zi — _L.U 


u .u 


fgas Gal2 


0.6-0.8 


expected at z=1.5 


0.8 


Orbit 


near polar 


to form the ring 




Gall incy b 


10 to 90 


Giant Stream 


35-75 


Gal2 incy b 


-30 to -110 


Giant Stream 


-55 to -110 


Gall incz c 


90 to 110 


Giant Stream 


90-110 


Gal2 incz c 


90 to 110 


Giant Stream 


90-110 


Spin Gall 


prograde 






Spin Gal2 


retrograde 


significant remnant disk d 




Tpericenter 


20-30 kpc 


see the text 


22-30 kpc 


Feedback 


high-medium 


to preserve gas 


high/ varying 


C-star 


0.004 f -0.03 


to preserve gas 


0.03 



"Daddi et al. 



12010) found / gas =0.5-0.65 in galaxies with M, 



baryon 



=0.8-2.2 IO^Mq at z=1.5 



"Orientation of the angular momentum of Gall relative to the orbital angular momentum, y axis 
c Orientation of the angular momentum of Gall relative to the orbital angular momentum, z axis 
d Following |Hopkins et al.||2008| > 

e In some simulations, the feedback is assumed to be high before fusion and later on, assumed to drop to the medium or low 

feedback values of |Cox et al.| ( p008^ ; see also section 4.2) 
•f Another way to preserve the gas before fusion in order to allow a significant amount of gas in the disk 



also adopted by |Cox et al. (2006); see also an ob- 
servational support from van der Kruit ( 2007) and 
references therein). We have to preserve enough 
gas immediately after the fusion to preserve the 
formation of a significant thin disk with more than 
50% of the baryonic mass. 

During our investigations to optimise the or- 
bital parameters, we run the GADGET2 code with 
150000 particles. In most simulations, we have 
ensured that the mass of dark matter particles is 
no more than twice that of initial baryon parti- 
cles, i.e., gas or stars in the progenitors. During 
star formation events controlled by the gas vol- 
ume density (and following the Kennicutt-Schmidt 
law), we have also limited the number of newly 
formed stars per gas particle to three. This is an 
important limitation to the simulation as we have 
verified that when newly born stars are relatively 
too small in mass, they are artificially scattered 
due to their encountering heavy dark matter par- 
ticles. Similarly, we need to keep the mass of dark 
matter particles low enough to avoid non-physical 



disruption of the newly reformed disk at the end 
of the simulation. In the following we have tested 
that the decomposition of the newly formed galaxy 
does not depend on the adopted number of parti- 
cles within the range of 150 000-800 000. 

However, simulating the first order properties 
of faint structures like the Giant Stream requires 
a high particle number: the Giant Stream is made 
of 3xlO 8 M , which corresponds to l/366th of the 
total baryonic mass of the galaxy, and lets only a 
few hundreds of particles in the Giant Stream for 
a 150 000 particles simulation. Having fixed the 
range of parameters to reproduce the observed de- 
composition of M31 in the bulge, and the thin and 
thick disks, we used a limited number of simula- 
tions with 300 000 to 800 000 particles to better 
identify the formation of faint structures. 

Table [l] lists the adopted range of parameters 
which have been investigated (2nd column) and 
adopted (column 4) for several models considered 
in this study. These models are all part of the same 
family of orbital parameters, and they differ essen- 
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tially by various adjustments of the star formation 
history or small variations of orbital parameters. 
We have begun with a model with moderate gas 
fractions (40% and 60% in Gall and Gal2, respec- 
tively) to verify whether a significant thin disk 
can be formed after fusion, assuming somewhat 
extreme conditions (e.g., high feedback and very 
low star formation efficiency) to prevent star for- 
mation before the fusion epoch. We then realised 
that most of the gas preserved before fusion forms 
a very thin disk (see also Abadi et al.|2003 ), which 
also hosts some of the pre-existing stars before fu- 
sion. However, the thin disk fraction is generally 
not much larger than the available gas fraction be- 
fore fusion. Because about 65% of M31 stars lie in 
the thin disk, we have generated for our current 
models a higher gas fraction that allows us to test 
less extreme conditions for feedback and star for- 
mation efficiency. Indeed, such gas fractions are 
quite common at z=1.4 (see Daddi et al.||2010" |, a 
redshift which corresponds to 9 Gyr ago, just be- 
fore the interaction between the M31 progenitors 
(see section 2) . Finally the models adopted in Ta- 
ble [I] (column 4) are developed to optimise the 
reproduction of the decomposition of the galaxy 
in the bulge, and the thin and thick disks, and the 
disk scale length (see sections 4.1 and 4.2), the 
presence of the 10 kpc ring (see section 4.3), the 
Giant Stream and its kinematics (see section 4.4), 
as well as the fraction of stellar ages in most of 
these substructures (see sections 4.2 and 4.4). An 
overview of the goodness of each model is provided 
in section 5. 

Figure [2] shows the different steps of the merger 
for one of our models. At the first passage (at 
T=0.7 Gyr), tidal tails are formed and enrich the 
halo, although they are mostly diluted at the fu- 
sion epoch. Part of the material ejected during the 
second passage and the fusion progressively return 
to the galaxy after the fusion. The material asso- 
ciated with the tidal tail pointing to the bottom 
of the galaxy (see panel at T=4.2 Gyr) possesses 
an angular momentum very different from that of 
the newly formed disk. This perturbation may 
explain such faint features like the Giant Stream, 
and could be persistent for several Gyr after fusion 
(see section 4.4). 



4. Results and comparison with observa- 
tions 

4.1. Decomposition of the newly formed 
galaxy in sub-components 

We choose to decompose the newly formed 
galaxy in three components, according to the rel- 
ative strength of the angular momentum in the 
direction of that of the thin disk, which is ap- 
proximately the orbital angular momentum. A 
very accurate determination of the thin disk an- 
gular momentum is possible because most of the 
preserved gas particles before the fusion naturally 



form a thin disk (see also Abadi et al. 20031 in 



which most of the star formation occurs after fu- 
sion. We have used the stars born 0.5 Gyr after 
the fusion to determine the angular momentum of 
the newly formed disk (see Fig. |3j. With this 
technique, we obtain an accuracy of ±1 degree for 
determining both the disk PA and inclination. 

Fig. [4] (top panels) shows the ratio of the an- 
gular momentum along the y axis to the total 
angular momentum {Jy I J total) as a function of 
the x-axis along which the thin disk is projected 
(see Fig. [3|. The fi nal disk rotates anti-clockwise 
(see |Chemin et ah 20091 and the thin disk (illus- 
trated in Fig. |3j) is easily recognisable in the Fig. 
[4] top-left panel for values ranging from -0.9 and 
-1, which corresponds to the gaseous thin disk. 
Other recently formed stars are located within a 
compact structure without preferred angular mo- 
mentum orientations, which corresponds to the 
bulge. Ilbata et al.l (120071) pointed out the con- 



fusion in the literature between the spheroid com- 
ponent and the thick disk. They revealed that the 
minor-axis region between projected radii of 7 kpc 
< R < 18 kpc is strongly affected by a rotationally- 
supported thick disk that corresponds to approx- 
imately one tenth of the thin disk stellar mas^J 
For a better comparison with observational data, 
we assume that the thick disk is made by all par- 
ticles (except those of the thin disk) whose angu- 
lar momentum is dominated by rotational motions 
along the thin disk. It means that Jy has to be 
larger than the combination of other components 



3 It is important to notice that this region (quoted as 
"spheroid and thick disk off-axis" in Fig. Q could well 
be dominated by the spheroid (Tom Brown, private com- 
munication), depending on the precise density profiles of 
both spheroid and thick disk. 
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of the angular momentum (i.e., \J J\ + J§ ) , i.e., 
Jy I J total should be larger than -l/y/2 =-0.707, 
as illustrated by the dot-dashed lines in Fig. [4] 
We have generated the decomposition proposed 
by Abadi et al. (2003) for comparison. In their 
scheme, the thick disk is the residue after subtract- 
ing the bulge and the thin disk assuming a bulge 
symmetrically distributed around Jy=0. While 
the thin disk fraction is found to be very simi- 
lar with both approaches, the residual thick disk 
from Abadi et al's method is almost twice that 
provided by our method. We notice that a signifi- 
cant fraction (up to 50%) of such a residual thick 
disk is made of particles which are not dominated 
by rotation around the Y-axis, which well explains 
the discrepancy. For comparison with the |Ibata et| 
al. ( 2005 1 observations of a rotationally-supported 



thick disk, we keep our decomposition as described 
above. 
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Fig. 3. — Distribution of newly formed stars 0.5 
Gyr after the fusion which occurs 5 Gyr after the 
beginning of the simulation, for a 2.8:1 merger 
with r per =30 kpc, Gall incy=65 and Gal2 incy=- 
90. Ten gigayears after the simulation, the diam- 
eter of the newly formed disk reaches ~ 36 kpc 
for a thickness of 4 kpc (see thin lines) in a model 
for which high feedback has been assumed during 
the entire duration of the merger. Here the thin 
disk has been projected along the X-axis, and its 
angular momentum is oriented along the Y-axis. 



Table 2: Decomposition in mass (unit=10 10 M Q ) of 
the newly formed galaxy 9 Gyr after the simula- 
tion, for various parameters of our modelling. For 
all models but the last one, high feedback has been 
maintained during all the simulation. For the last 
simulation, feedback has been abruptly dropped 



from high to low values (from Cox et al. 20081, 
3.5 Gyr after the beginning of the simulation and 
~ 0.7 Gyr after the fusion. 



Parameters 


Comp. 


Thin 


Bulge 


Thick 






disk 




disk 


Tpericentei ^4 


Stars 


2.66 


2.51 


0.85 


mass ratio =3.5 


Gas 


3.24 


0.15 


0.09 


Gal inc=65 & -110 


Fraction 


62% 


28% 


10% 


Tpericente'i 24.8 


Stars 


2.72 


2.31 


0.78 


mass ratio =3.0 


Gas 


3.34 


0.12 


0.11 


Gal inc=65 & -89 


Fraction 


65% 


26% 


9% 


fperic&ntei 24.8 


Stars 


2.38 


2.67 


0.95 


mass ratio =2.8 


Gas 


2.95 


0.13 


0.10 


Gal inc=65 & -89 


Fraction 


60% 


29% 


9% 


Tpericentei ^4 


Stars 


2.25 


2.84 


0.88 


mass ratio =2.5 


Gas 


2.91 


0.12 


0.09 


Gal inc=65 & -89 


Fraction 


57% 


32% 


10% 


Same as above 


Stars 


3.62 


2.61 


1.22 


feedback 


Gas 


1.32 


0.07 


0.12 


changed at 3. 5 Gyr 


Fraction 


55% 


29% 


14% 



4.2. Thin and thick disks: feedback pre- 
scriptions 

Table [2] describes how the decomposition de- 
pends on the choice of orbital parameters. From 
the 2nd to the 4th sets of parameters, only the 
mass ratio varies and the bulge fraction increases 
as the mass ratio decreases as shown by[Hopkins ct 
al. (2009). However, we have also tested a higher 



mass ratio (1st set of parameters) with a differ- 
ent set of inclinations for the progenitors. Thus 
the bulge fraction may also depend on other or- 
bital parameters, here the relative inclination of 
the progenitors. For example, it reaches a max- 
imum when the difference of inclination between 
the two progenitors is close to 180 degrees, possi- 
bly because more stars take radial orbits and fall 
into the bulge. It will be shown later in section 4.4 
that this resonance has a more considerable effect 
on the matter that is stripped from the satellite to 
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Fig. 4. — Top: adopted decomposition of the 
newly formed galaxy (same model as in Fig. |3j, 
10 Gyr after the simulation. Here the thin disk 
has been projected along the X-axis as in Fig. [3] 
Blue, green and red dots shown in the three pan- 
els represent stars formed 0.5 Gyr after the fusion, 
stars formed between the first passage (0.5 Gyr 
after the beginning of the simulation) to 0.5 Gyr 
after the fusion, and stars within the progenitors, 
respectively. The two dot-dashed lines delineate 
the thin and the thick disks, from the bottom to 
the top. Bottom: particles distribution in the thick 
disk (left panel) and in the bulge (right panel). 
Here the thick disk and the bulge have been pro- 
jected at the right orientation of M31 thin disk 
(e.g., PA=38 degree and inclination=77 degrees). 



the tidal tail at the 2nd passage. 

Almost all the simulations we have generated 
are re-forming significant thin disks containing 
more than half of the baryonic mass. This is due 
to our feedback prescriptions that always preserve 
a significant fraction of gas just before the fusion. 



As stated by Cox et al.| ( |2008[ ), very little is known 
about the requisite conditions enabling the star 
formation to occur. An important prerequisite, 
however, is the presence of stable gas rich galactic 
disks at high redshifts, that could be the progen- 
itors of M31. |Cox et all (12008} have tested feed- 



back for an isolated disk galaxy with an initial gas 
fraction of / 9as =0.52. They found that the gas 
fraction is mostly preserved within one Gyr if me- 



dian to high feedback prescriptions are adopted, 
while low feedback decreases the gas fraction to 
0.37 in the same amount of time. 

In our current models we have assumed f gas ~ 
0.65: the re-formation of a significant thin disk 
implies that most of the gas is not transformed 
into stars before fusion. This calls for relatively 
high feedback, at least intermediate between the 



median and high values of Cox et aL] ( |2008[ ) We 
have verified that for half the value of high feed- 



back of Cox et al. ( 2008 ) , most of the initial gas is 
preserved before the fusion. However, re-forming 
a significant stellar thin disk requires a significant 
change of the feedback at or after the fusion, in 
order to transform stars from the gas in the newly 
reformed disk. It results that from the merger of 
realistically gas- rich galaxies at 1.4, we may 
reform a galaxy resembling M31, with the condi- 
tion of a transitory history of the feedback dur- 
ing the merger, from high to low values. There 
could be some theoretical grounds favouring such 
a history, however. For example, in a gaseous rich 
and almost pristine medium, first supcrnovae are 
likely massive and generate high feedback as they 
may delay star formation for up to 100 Myr (e.g., 
Bromm et al. 2009 1 . The fusion of two gaseous- 



rich galaxies corresponds to a severe mixing of 
most of their components, accompanied by large 
star formation rates. Then the metal abundance 
in the remnant becomes larger and more homoge- 
neously distributed than it was in the progenitors, 
before fusion. The occurrence of very massive, 
primordial supernova is unlikely in such a mixed 
medium leading to a possible transition from high 
to low feedback. Before the fusion there are many 
fewer exchanges between the two interlopers, and 
it might be realistic to assume a negligible change 
in the feedback history. 

Of course, the above is mostly made of con- 
jectures in the absence of observations of gas-rich 
galaxies at high redshift showing no or very small 
amounts of star formation, i.e., the high-redshift 
counterparts of the present-day low surface bright- 
ness galaxies. We may, however, verify whether 
our adopted star formation history is consistent 
or not with the star formation history revealed in 
each M31 substructure displayed in Fig. [T] The 
thick disk is indeed mostly made by matter return- 
ing from the tidal tails. It is comprised of approx- 
imately 80% of stars older than 8 Gyr, the age of 
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other stars being mostly from 5 to 8 Gyr. Such a 
distribution is well reproduced by our model with 
high feedback values (86% of old stars) or with half 
this value (80%). Because the thick disk shows 
a similar fraction of stars in the two bins set by 



Brown et al. (2008) at 5.5 and 7 Gyr, it is proba- 
ble that there are no significant changes of the star 
formation history between the two corresponding 
epochs, i.e., between the first and second passages. 

Within these prescriptions our models predict 
quite well the fraction of each sub-component of 
M31 (see Table [2]). Note also that the modelled 
thin disks show rotation curves and scale lengths 
that are in good agreement with observations (see 
Fig. [5]). In the following sections we examine 
whether this model could also reproduce other 
large-scale structures such as the ring and the Gi- 
ant Stream. 

4.3. The 10 kpc ring 

The decomposition of the light profile was done 
by fitting the radial profile with a 3 components 
model (bulge, thin, thick disks). This provides 
a thin disk with a Sersic index of 1 and a scale- 
length varying from 4.8 to 5.1 kpc for models with 
mass ratio ranging from 3 to 2.8. Such values are 
slightly smaller than the 5.8±0.4 kpc that was 
adopted by |Hammer et al.| ( |2007[ ) for the M31 
disk. All these decompositions (see Fig. [5]) re- 
veal the presence of a prominent ring at 10 kpc 
which corresponds well to the observed ring. Sim- 
ilarly all models reproduce the 10 kpc gaseous ring 
(see Fig. [6]) as seen in HI observations. This is 
not unexpected as all the orbits are close to po- 
lar which favours the formation of such prominent 
and persistent structures. In fact, HI observations 
do not provide many constraints on our modelling. 
They only discard parameters for which the disk 
is too warped, i.e., implying that the difference 
between the two inclination angles of the progeni- 
tors should not be too close to 180 degrees. More- 
over, we notice (see Fig. [6| that with a constant 
high feedback during the simulation, there is still 
gas within the bulge of the remnant galaxy. How- 
ever this gas is consumed in models that assume 
a transition to low feedback immediately after the 
fusion. 
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Fig. 5. — 3:1 merger with r per j cen t er =24.8 kpc, 9 
Gyr after the beginning of the simulation. Top: 
the rotation curve of the thin disk compared to 
Carignan et al. (2006). There is a particularly 
good agreement between our modelling and data 
from Carignan et al. (2006). Bottom: the decom- 



position of the mass profile in 3 components which 
evidences the lOkpc ring and provides a thin disk 
scalelength of 5 kpc. 

4.4. The Giant Stream and its kinematics 

The Giant Stream is a very faint stellar struc- 
ture, which extends from 50 to 80-100 kpc from 
the M31 centre. We have modelled it as being 
caused by particles coming back from the tidal 
tail formed just before fusion. Fig. [2] describes 
the formation of such a structure that is aligned 
along the trajectory of the satellite which falls into 
the mass centre at the fusion, 4.5 Gyr after the be- 
ginning of the simulation. We verified that, within 
the family of orbits we choose, the strength of the 
tidal tail (see Fig. [2] panel at 4.2 Gyr) depends 
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properties to identify peculiar structures able to 
persist after the remnant phase (see Fig. [7|. 





Fig. 6. — Comparison of the HI gas from Braun 



et al. (20091 with modelling of gas particles for 
3:1 mergers at T=9 Gyr. On the top-right and 
bottom-left models are with r per =24.8kpc, and inc 
of 65 and -95 degrees for the main galaxy and the 
satellite. The two models have a feedback vary- 
ing from high to low values near the fusion time, 
and they do not show gas excess in the central 
bulge. They only differ by a small variation of 5 
degrees in their inclinations along the z-axis. On 
the bottom-left is shown a model with a constant 
high feedback and with a difference between the 
two galaxies inclination of 170 degrees. It illus- 
trates that changing feedback is necessary to ex- 
plain the absence of gas at the centre and that 
inclinations should not be too close to 180 degrees 
to avoid too large distortions of the final gaseous 
disk. 

on the inclination of the progenitors relative to the 
orbital angular momentum (see Table]!]). Optimal 
values are found for large values of the inclination 
of the main galaxy (> 50 degrees) and especially 
for differences between the inclination of the two 
progenitors between 140 to 170 degrees, i.e., not 
too far from the resonance at 180 degree. 

Modelling the Giant Stream is not an easy task 
as the structure is very faint compared to other 
large substructures of M31, even the thick disk. As 
in section 4.1 we make use of angular momentum 





. . 1 , 1 . . . , , , 1 , 1 , . . 

INTERMEDIATE & - 




'. OLD STARE 


I 






1 1 



50 TO T70kpe 



dijat tali of r 
all particles main (AM 
near disk) 





— m — 

RETURNING 
PARTICLES FROM 
THE TIDAL TAIL 
OF THE 
SATELLITE 



100 



LOO 
, Wee) 



100 200 300 



Fig. 7. — 2.8:1 merger with r pericenter =30 kpc. 
The distribution of the angular momentum 5 Gyr 
after fusion using spherical coordinates, assuming 
the thin disk perpendicular to the y-axis. Top-left: 
stars formed before fusion (intermediate and old 
stars) which show a broader distribution than thin 
disk stars as it is expected for the thick disk that 
hosts most of these stars (see Table [2J. Top-right: 
angular momentum distribution for particles be- 
yond 50 kpc from the centre of the remnant. Most 
particles in the outskirts have angular momentum 
following the two tidal tails, one associated with 
the main progenitor, the other associated with the 
satellite. Bottom-left: same after removing both 
thin and thick disks: only the particles in the loops 
are kept. Bottom-right: the distribution of angu- 
lar momentum for all the particles associated with 
the tidal tail of the satellite, evidencing that these 
particles are lying in the $=225 degrees plane. 

Fig. [TJshows the distribution of the angular mo- 
mentum for the thick disk (top-left panel) and at 
the galaxy outskirts, 5 Gyr after fusion. The top- 
right panel shows two peaks which correspond to 
the two tidal tails formed after the second passage. 
The tidal tail associated to the main galaxy (see 
Fig. [2] at 4.2 or at 4.5 Gyr) is spiralling around 
the thick disk and its angular momentum peaks at 
<f>=290 degrees, which is indeed part of the thick 
disk (see top- left panel). This peak disappears 
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in the bottom-left panel for which the thin and 
thick disks have been removed. The residual par- 
ticles form another well identified peak at $=225 
degrees. We have verified that they correspond 
to particles returning from the tidal tail of the 
satellite: in the bottom-right panel, selected par- 
ticles of the tidal tail of the satellite indeed show a 
very similar distribution with a peak at $=225 de- 
grees, which is offset by 45 degrees from the disk 
angular momentum. The selection of these par- 
ticles has been done in a rather crude way (see 
Fig. [2| panel at T=4.2 Gyr): we simply "cut the 
tidal tail" at the time of its formation, up to the 
edge of the satellite. However this crude selection 
misses some particles, for example those which are 
stripped from the satellite at later phases (e.g. 3rd 
passage, see Fig. [2] panel at 4.5 Gyr). Comparing 
the counts of particles having an angular momen- 
tum around the $=225 degrees peak, we find that 
our selection recovers 90%, 70%, 40% and 15% 
within shells of 100-170kpc, 50-100kpc, 30-50kpc 
and 20-30kpc, respectively. As expected, particles 
at the very edge of the tidal tail are the most diffi- 
cult to pre-select. Five gigayears after fusion those 
particles have returned to trajectories close to the 
galaxy centre. 

Tracing the trajectories of selected particles in 
the tidal tail is very instructive: they correspond 
to an important component of the angular mo- 
mentum which is at odd with the rotation of both 
the thick and thin disks (see Fig. [8|. 

Fig. [H] describes the temporal evolution of the 
trajectories of tidal tail particles returning to the 
galaxy. When falling to the galaxy centre, the par- 
ticles are looping around the galaxy with a peri- 
center at ~ 25 kpc. These particles are describ- 
ing several elliptical orbits -or loops- around the 
galaxy potential. We identified 4 of these loops, 
although they might be more numerous. Fig. [8] 
shows that with increasing time, they show a pre- 
cession around the galaxy as well as they are grow- 
ing. This is because stars coming back from the 
tidal tail are returning from higher elevation at a 
later epoch and are coming back with increasing 
energy. We also note that the loops are somewhat 
thick and, as such, they could better be called 
tores. 

The permanent rain of stars from the tidal tail 
ensures very long-term streams of stars within the 
galaxy outskirts well after the remnant phase. In- 



falling stars are trapped into elliptical orbits with 
the remnant galaxy right at the foci, which make 
the orbit stable. These orbits are those expected 
for extremely small satellites (single stars!), in the 
absence of tidally induced forces. Loops are ex- 
pected to persists for several billion year^] after 
fusion and, moreover, they are permanently fed 
by new particles coming from the tidal tail. These 
streams and loops are all within a common plane, 
45 degrees from the PA of the rejuvenated disk as 
it is evidenced by Fig. [8] (panels h and i) . This is 
our proposed model for the formation of the Giant 
Stream which indeed points toward the centre of 
the galaxy. Let us now examine how this model is 
compliant with several observations of the Giant 
Stream properties. 
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Fig. 9. — Left: same model as in Fig 3. The 
distribution of heliocentric velocities in the Giant 
Stream region at distances indicated in each panels 



with references to fields from Gilbert et al. ( 2009 1 . 
The dot-dashed line reproduces the black line of 
Figure [l] of Ibata et al. ( 2004 1 representing the 
high negative velocity edge. Green dots represent 
the mass weighted average for each panel. Right: 
same for particles returning from the tidal tail. In 
this plot we have accounted for several snapshots 
in the simulation to artificially increase the num- 
ber of stars in the central regions (see the text). 

The kinematics of stars follows the trend evi- 



denced by Ibata et al. (2004) and Gilbert et al 



( 2009 ) : stellar particles in the streams reach large 



negative values of their heliocentric velocities to- 

4 We have verified that the resulting dark matter density 
profile in the remnant is slightly oblate (e=0.2) shortly af- 
ter the fusion of the two galaxies and then stable: then it 
could not strongly affect the trajectories of returning star 
particles from the tidal tail. 
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Fig. 8. — 2.8:1 merger with r per i cente r—SO kpc. From panels (a) to (g): Snapshots in the (y, z) plane of 
trajectories for particles (red dots) which are part of the tidal tail of the satellite. Panel (a) shows the 
system just before fusion when particles are stripped from the satellite. Panels (b) to (f) show the particle 
motions for time intervals of 0-lGyr, l-2Gyr, 2-3Gyr, 3-4Gyr and 4-5Gyr after fusion, respectively. Panel 
(g) shows the distribution of particles 5.5 Gyr after fusion. For illustration, we have added in all panels the 
thin disk (black) as it is, 5.5Gyr after fusion. It shows the formation of the different loops which are drawn 
by particles when coming back from the tidal tail. Panels (h) and (i): Same as panel (g) but projected in 
the observer's frame in the (x, y) directions illustrating that the loops are inserted within a thick plane at 45 
degree from the disk PA. Panel (h) shows the same particles as panel (i) for which the heliocentric velocities 
are coded following [Chemin et"al] ( |2009| ), i.e. cyan (<-500 km/s), blue (-500 to -400 km/s), green (-400 to 
-200km/s), red (-200 to -100 km/s) and magenta (> -lOOkm/s). It illustrates a disk rotation very similar 
to the observed one as well as the velocity of the Giant Stream which is reaching larger negative velocities 
when it reaches the M31 centre. 



wards the centre (see panel i of Fig. [8]). Fig. [9] 
shows the distribution of heliocentric velocities for 
stellar particles selected in the region of the Giant 



Stream within & b s from -90 to -140 degrees. Ob- 
servations (see_^ig 1 Jl]of Ibata et al. (2004) and 
Fig. [5] of |Gilbert et al.| ( |2009| )) show a broaden- 
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Fig. 10. — Comparison of observations (top-left 
from Ibata et al.||2005| with our modelling; at 
top-right, we adopted a mass ratio of 2.8:1 with 
r per —30 kpc observed after 10 Gyr, at bottom-left, 
a mass ratio of 3:1 with r per —25 kpc observed af- 
ter 9 Gyr. The bottom-right panel shows the same 
as that of the top-right, except that we have in- 
verted all axes. 4 to 10 snapshots have been used 
to optimise the contrast. 



ing of the velocity distribution towards the galaxy 
centre that is well reproduced by our simulations. 
Besides this, our model also well reproduces the 



distribution of stellar ages observed by Brown et 



al] (|2008), as we find 80%, 20% and 0% of stars 
older than 8 Gyr, from 5 to 8 Gyr and younger 
than 5 Gyr, respectively. These fractions can be 
regulated by changing the feedback prescriptions 
(see section 4.2). 

We have attempted to reproduce the photomet- 
ric properties of the Giant Stream. The difficulty 
of such an attempt is due to the fact that the Gi- 
ant Stream is very faint when compared to other 
sub-structures of M31, which is a problem also en- 
countered by observers. Even with 500 000 parti- 
cles, the number of stellar particles in the Giant 



Stream is below 1000 which obviously limits a de- 
tailed reproduction of the structure. Within an 
angle of Q b s from -90 to -140 degrees and from 25 
to 80 kpc, the total mass of the Giant Stream is 4.8 
and 4.4 10 8 Mq in the two examples illustrated in 
Fig. 10 Using other variations of the model pa- 
rameters we find a range from 2 to 5 x 10 8 Mq , 
which easily includes the observational value of 3 
10 s M Q (see |Mori fc Rich||2008[ ). 

All our models (column 4 of Table [T]) show 
streams of stars which belong to a plane that in- 
tersects the observational plane in a direction 45 
degrees from the thin disk PA. We verify that to 
optimise the reproduction of the observed features 
in position and strength, we can use two differ- 
ent kinds of symmetries. First, we can rotate the 
whole system around an axis perpendicular to the 
thin disk, as this would not affect the galaxy de- 
composition in substructures, nor the kinematics 
of thin and thick disks. Second, we hnd that grav- 
ity is not sensitive to a complete inversion of the 
initial (x, y, z) coordinate system. This has the 
advantage of being able to reproduce the Giant 
Stream with many variations of our model param- 
eters, but it also has the inconvenience of provid- 
ing a very large amount of space parameters to 



investigate. Fig. 10 illustrates some examples of 
this exercise for two models with 500K particles 
for which the Giant Stream mass is slightly larger 
than the observed one. Much larger number of 
particles are certainly required to reproduce at the 
same time all the morphological details of the thick 
disk, although we generally reproduce quite well 
the NE and W shelfs. 

5. Discussion and Conclusion 

By construction, all models reproduce quite 
well the distribution of stellar ages in the vari- 
ous substructures, as they are described in Fig. [I] 
This is especially true for both the thick disk and 
the Giant Stream (see sections 4.2 and 4.4) for 
which our predicted compositions of stellar ages 
are from 10% to 20% of intermediate age stars 
(5-8 Gyr) , the rest essentially being made of older 
stars. However, we experienced some difficulties to 
reproduce the distribution of the 21kpc and 35kpc 
"Halo" fields oflBrown et all p)08|). We indeed 



find a fraction of intermediate age stars ranging 
from 5% to 15% , while according to Figure [TJ it 
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should be near or below 5%. It would be tempting 
to use this property for selecting the best model 
of M31. However, in our models intermediate age 
stars are related to low mass particles, three of 
them being produced by a single gas particle. It 
may lead to an artificially scattering of those in- 
termediate age star particles (see section 3) , which 
could be only solved by simulations with a much 
larger number of particles than 500K. 

All models as defined in Table [TJ column 4 can 
reproduce with a quite good accuracy the bulge, 
and thin and thick disk fractions, as well as the 
ring and the shape, the mass, and the kinematics 
of the Giant Stream. The examples shown in Fig. 
[TOl are cases for which the Giant Stream stellar 
content is slightly larger than the observed one, for 
illustration purposes. We also note that the shapes 
of the predicted thick disk and Giant Stream can 
significantly change -within a given set of orbital 
parameters (see Table [I])- either because they can 
be modified through a rotation around the disk 
axis or because they could fluctuate with the as- 
sumed time after the beginning of the simulation. 
It is beyond the scope of this paper to search for 
the best model that reproduces all the details of 
M31 substructures. This would require launching 
several hundreds of simulations with several mil- 
lions of particles. We have, however, learned from 
this analysis that: 

• A thin disk such as that of M31 can be repro- 
duced with a combination of ~ 65% of gas 
in the progenitors and a star formation his- 
tory that prevents gas consumption before 
the fusion; low efficiency of star formation is 
thus predicted in primordial gaseous phases, 
while when the medium has been sufficiently 
enriched and mixed during fusion, the gas 
is more easily consumed; this transition has 
essentially been modulated through the as- 
sumed feedback history; 

• The same star formation history also re- 
produces a rotationally-supported thick disk 
with a stellar composition similar to what is 
observed; and 

• also fractions of bulge/thin disk and thick 
disk similar to what is observed within 10% 
accuracy. 



• These polar orbits always reproduce the ob- 
served lOkpc ring, and in most cases, the 
gas map is very similar to the observed one 
(see Fig. [6]), especially when the gas in the 
galaxy centre is allowed to form stars after 
the galaxy fusion; 

• The predicted Giant Stream is made of stel- 
lar streams due to returning stars from 
a tidal tail formed prior to fusion from 
stripped particles of the satellite; these stars 
are drawing streams along loops which are 
inserted into a thick plane that is 45 degrees 
inclined relative to the disk PA (see Fig. [8]). 
Such streams are very long-lived because re- 
turning particles are trapped into such loops 
for periods larger than several gigayears, and 
also because the streams are permanently 
fed by new particles falling back from the 
tidal tail at all epochs. 

The success of modeling M31 properties as re- 
sulting from a past major merger event is patent. 
Its predictive power is however limited by the com- 
plexity of the system as well as the very demand- 
ing number of particles required to reproduce in 
details all substructures of M31. An important 
limitation could also be the accuracy of the calcu- 
lations, which requires tracking with an unprece- 
dented accuracy positions and velocities of parti- 
cles for almost a Hubble time. In fact, we cannot 
claim to have determined the best major merger 
model of M31, although it is likely within the lim- 
its discussed below. 

We did not either try to investigate the whole 
parameter space to which the above results hold, 
although examination of Table [T] led us to sus- 
pect that it is likely large. However some specific 
links between parameters and observations may 
provide some empirical limits. The bulge frac- 
tion can be reached within a 10% accuracy for 
all models with the mass ratio of 3±0.5. In the 
frame of our model for M31, the halo is enriched 
by stars coming from the progenitors before the in- 
teraction, while the thick disk and Giant Stream 
are further enriched by stars formed before the 
first passage and fusion. This (see Figure [I] and 
section 2) is compliant with pericenter radii of 
25±5 kpc. Further modelling and measurements 
of stellar ages in several fields in the outskirts of 
M31 would certainly help to provide more accu- 
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rate values. This would also provide better accu- 
racy for the important epochs of the interaction 
which would have begun 8.75±0.75 Gyr ago with 
a fusion time 5.5±0.5 Gyr ago. There is still work 
needed to properly constraint other orbital param- 
eters such as the inclinations of progenitors rela- 
tive to the angular momentum direction. We are 
far from having investigated the whole range of or- 
bital properties which produces a strong tidal tail 
from stripped material from the satellite, and then 
a structure similar to the Giant Stream. Table [3] 
describes all the models presented in this paper, 
which are part of about 50 tested models. Fig. [6] 
demonstrates that models with differences of pro- 
genitor inclinations near 180 degrees show not only 
a maximum of resonnance (thus a very strong tidal 
tail) but also galactic disks which are probably too 
warped or distorted. We have mostly investigated 
models with the main galaxy inclination of 55 to 
70 degrees, but there could be other solutions for 
different angles. All the models in Table |3j as well 
as models in the range described by Table [l] (4th 
column) reproduce quite well the observations, in- 
cluding the Giant Stream. A much more detailed 
analysis is necessary to determine the best model 
parameters, and this could be done by comparing 
in detail the constraints provided by the gaseous 
disk, the thick disk, and the Giant Stream shapes 
as well as by any other structures discovered in 
the M31 halo field. 

While writing this paper, we become aware of 
the discovery of many other structures in the halo 



of M31 (Mackey et al. 20101, especially after the 



successful deep imagery of the whole field sur- 



rounding M31 by the PANDA team ( McConnachie 



et al. 2010). This magnificent image reveals the 



presence of loops surrounding M31 which strik- 
ingly resemble our model predictions (see Fig. 
[8]). In particular the northern loop found by the 
PANDA team which extents up to 120 kpc from 
the galaxy centre provides particularly strong con- 
straints on our model. Indeed the first loop in 
Fig. [8] can accurately reproduce the northern 
loop (called NW stream in the incomplete image 
of Mackey et al. 20101, and other loops are well 
matched with the Eastern arc and Stream C. Fig. 
[TT] shows two examples of models that succeed in 
reproducing most of these features along with the 
Giant Stream and most other substructures (thin 
disk, bulge, and thick disk) of M31. In fact we also 
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Fig. 11. — Two models for which most features 
have been reproduced including those from the 
new PANDA survey. Top: a model with 500K 
particles assuming a 3:1 merger, r per =25 kpc, and 
inversion of the initial (x, y, z) coordinate sys- 
tem. On the left are shown the returning particles 
from the tidal tail as in Fig. [8j panels (h) and 
(i), with black dots and with the rotating disk in- 
serted (same code as that of Fig. |8| . On the right 
all particles of the simulation are shown, includ- 
ing 4 snapshots around T=8.72 Gyr. This model 
is with Galli ncy =65, Gal2 incy =-95, Gall incz =85, 
Gal2i ncz =85, half the maximal feedback until 
0.5Gyr after fusion, then low feedback. Bot- 
tom: a model with 1M particles assuming a 2.8:1 
merger, r per =30 kpc and inversion of the initial 
(x, y, z) coordinate system. Here only 2 snapshots 
around T=9.05Gyr are shown in the left image. 
Here a constant high feedback history has been 
assumed and other parameters are Galli ncy =65, 
Gal2 incy =-89, Gall incz =U0, Gal2 incz =110. The 
Giant Stream mass includes 2.5 to 2.8 10 s M@. 



verified that the Giant Stream kinematics is still 
reproduced. Fitting these additional structures al- 
lows us to fix the rotation angle around the thin 
disk axis as well as solving the initial conditions of 
the system because only the inversion of the initial 
(x, y, z) coordinate system is able to reproduce the 
ensemble of features in Fig. [TT] The size of the 
loop also depends on the time elapsed after fusion 
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Table 3: Summary of model parameters from Table [T] used throughout the paper (all within the range of 
parameters of Table [Tj column 4). Angle, pericenter radii, number and mass of particles are given in degrees, 
kpc, 1000s of particles and (1O 6 M0), respectively. Feedback adopted values are (1): high feedback; (2): 5 
times medium feedback; (3) high feedback before fusion then low feedback and (4): 5 times medium feedback 
before fusion then low feedback. Model parameters for Fig. |4j [7J [HJ [9] and [10] are the same than for Fig. [3] 
the other model used in Fig. [TU] is the same as that of Fig. [TT^ i. The 6 last lines of the Table provides a 
check of the goodness of each model within ±20% of the observed value: the decomposition of M31 in the 
bulge, and the thin and thick disks (e.g. B/T=0.28±0.05), the disk scale length (e.g. 5.8±1.1 kpc), the 10 
kpc ring, the Giant Stream (e.g. 3±0.6 10 8 Mq), the stellar ages in the thick disk and in the Giant Stream 
(assuming fraction of > 8 Gyr stars to be 85±8%) and gas fraction in the galaxy within 30 kpc (7±2%). 
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"Models with a constant star formation with low efficiency (or high feedback) could not transform enough gas into stars after 
the fusion leading to a more gas-rich disk than observed; assuming a varying star formation (more efficient after the fusion) 
suffices to correct the discrepancy. 



and the beginning of the interaction is expected to 
be 8.75±0.35 Gyr within the range of parameters 
in Table □ 

As a whole, M31 and the complex structures in 
its halo can be reproduced by assuming a single 
major merger which may have begun 8.7 Gyr ago. 
The advantages of such a solution for M31 com- 
pared to the hypothesis of numerous minor merg- 
ers are as follows: 

• It explains most of the complexities of M31 
by a single event rather than numerous 
and minor merger events which have to be 
adapted in a somewhat ad-hoc way to each 
feature discovered in the M31 halo, including 
the thick disk; 

• It naturally explains stellar ages and metal- 
licities of most substructures in the M31 halo 
as well as why they show so many similari- 



ties; 

• It overcomes the increasing difficulty of iden- 
tifying the residuals of satellites assumed to 
be responsible for the different features in 
the M31 halo; 

• It is consistent with the kinematics of the 
M31 globular system (e.g., |Bekki~||2010[ ). 



Re-formation of disks is expected after gas-rich 



major mergers (Barnes, J.E. 2002 Springel et al 



2005b). It is more than just a curiosity since the 



disk rebuilding scenario has been proposed as a 



channel for the general formation of spirals ( Ham 
mer et aL]|2005 ), including specifically the case of 
M31 ( Hammer et al.| [2007). Gas-richness of pro- 
genitors is certainly high, as it has been estimated 
for many high-z galaxies. The requirement of a 
star formation with a low efficiency before the fu- 
sion can be reached through any mechanism that 
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Fig. 12. — Fe/H distribution of progenitor stars 
in the Giant Stream (red dashed line) and in the 
northern loop (blue solid line), for the same model 
as that of Fig. [TTfc ,. The absence of metal rich stars 
(e.g., -0.4<Fe/H<-0) is due to the fact that parti- 
cles returning from the tidal tail are coming from 
the outskirts of the progenitors. The absence of 
stars richer than Fe/H=-0.7 could explain why the 
northern loop is not detected with such a metallic- 
ity selection, while it becomes prominent for metal 
poor stars (Fe/H<-1.4). 



prevents star formation in relatively primordial 
medium. It could be caused either by higher feed- 
back efficiency in primordial medium as assumed 
in section 4.2, or because the gas in the progeni- 
tors is less concentrated than in present-day spi- 
rals as it is in present-day low surface brightness 
galaxies, or because cooling is less efficient in a 
relatively pristine medium, or a combination of 
all these factors. More generally, the expected in- 
crease of the gas metal abundance (expected to 
be slow before the fusion but very efficient dur- 
ing the fusion) may help to increase the molecu- 
lar gas fraction, the optical depth of the gas and 
the radiation pressure effects, all contributing to a 
change in the star formation history during the in- 
teraction (TJ Cox, private communication). The 
above is quite well illustrated a contrario by |Bour^| 
naud et al. (2010) who assumed that progenitors 



of high-z major mergers are very strong star form- 
ing, turbulent, gas-rich disks: they naturally find 



that remnants cannot be disk dominated within 
such conditions. 

There is still a considerable need to improve 
the modelling of M31, by using a much larger 
number of particles as well as by providing more 
kinematical and stellar age constraints in the nu- 
merous fields of view surrounding M31. We also 
need to reproduce with more accuracy the shape 
of the Giant Stream and of the thick disk con- 
tours. Metal abundance enrichment can also be 
implemented to more accurately fit the observa- 
tions. For example, the striking differences be- 
tween features at different metallicity in the halo 
found by the PANDA team cannot be retrieved 
by our modelling in the absence of such an imple- 
mentation^] We have attempted to verify whether 
such a trend is predicted by our model. At first or- 
der, one may "paint" each star in the progenitors 
assuming a reasonable metallicity gradient . We 



adopt the approach by Rocha et al. ( 2008 ) assum 



ing that metallicity gradient follows the dust gra- 
dient, with a scalelength 1.4 times larger than that 
of the stars. The range of metal abundance (Fe/H) 
is adopted from to -2.5. The stellar metal abun- 
dance of the satellite is likely significantly smaller 
than that of the main progenitor, because their 
stellar mass ratio (for our 3:1 models, see Table [TJ 
is 5. Little is known about the shape of the metal- 
stellar mass relationship at z ~ 1.5, although a 
serious attempt has been made at z^ 0.7 by |Ro-| 
drigues et al. (2008). Assuming this holds up to 



z=1.5, stars in the main progenitor are on aver- 
age 0.4 dex more metal rich than those from the 
satellite, a value smaller than the large scatter of 
the Rodrigues et al. (20081 relationship. Figure 



[12] shows the metallicity distribution of "painted" 
stars in both the Giant Stream (defined as in sec- 
tion 4.4, from 25 to 80 kpc and $ (, s from -90 to 
-140 degrees) and in the northern loop (defined 
from 55 to 130 kpc and & b s from to 90 degrees) 
discovered by the PANDA team. Because most 
of the stars in the northern loop are coming from 
the satellite, this could explain why this feature is 



5 This may also help to better reproduce the 35kpc halo field 
of |Brown et al.| l2007|l which is in-between the Eastern arc 
and the stream C jMackey et ah] 



■ |2010[ ): this also points out 
the need for resolving stellar streams and then for a better 
accuracy in following positions and velocities of particles; 
a similar need is for identifying carefully the particles in 
loops such as that proposed for the Northern loop. 
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not detected in the metal-rich map of the PANDA 
team (see Fig. 12 1. Such a prediction may stand 



unless intermediate age stars (those formed before 
the fusion) strongly affect the metal abundance of 
both the Giant Stream and the northern loop. On 
the other hand the same model correctly predicts 
the (observed) small fraction of such stars in the 
Giant Stream. 

Another important improvement in this mod- 
elling would be to consider the possible interaction 
of M33 which has been proposed to be responsi- 
ble for several features in the SE of the M31 halo 



by McConnachie et al. (2009). Although the dark 



matter density profile does not significantly affect 
our results (results shown in Fig. |j |3| [§ [5] [7) [sj [9) 

and 11 are similar when using a Hernquist profile 
instead of a core profile), we also need to verify 
whether these are affected by adopting different 
baryonic mass fractions, especially for lower val- 
ues than our adopted 20%. Indeed, a considerable 
change of the dark matter fraction may prevent 
stars from being stripped by the collision and then 
the distribution of baryonic matter in the outskirts 
may be significantly affected. 

However, we believe the result of such a model 
is very encouraging because it fulfills the Occam's 
razor principle: a single event may explain most -if 
not all- the properties of M31 and of its outskirts. 
It also proposes a mechanism for the Giant Stream 
which is simply fed by stars which are captured 
into loops orbiting around M31 for elapsed times 
which may reach a Hubble time. If M31 is actually 
the result of an ancient major merger, there will be 
a considerable need to revise our knowledge about 
our immediate environment and in cosmolog}j^] for 
the following reasons: 

• Merging of gas-rich distant galaxies can eas- 
ily produce large thin disks by assuming a 
less efficient star formation before the fusion; 
this simple scheme supports the rebuilding 
disk scenario for many giant spiral galaxies 
at least for those sharing similar properties 
(enriched halo stars) as M31; 



3 The history of galaxy formation theory has been strongly 
influenced by our knowledge of the Milky Way and the the- 
ory of spiral galaxy formation a la Eggcn ct al. ( 1962 ) is still 
lively; the fact that the second closest galaxy, M31, could 
plausibly be a major merger may remind us that the theory 
of galaxy formation could require significant adjustments. 



Most models show the formation of tidal 
dwarves, including those from tidal tails 
formed after the first passage, and some of 
them could be part of the satellite system of 
M31; their properties, after comparison with 
observations, may be used in the future as 
a constraint to the orbital parameters in 
further modelling; 

Up to 15% of the material is ejected after 
the fusion and this may populate the whole 
Local Group including in the direction of 
the Milky Way; this leads Yang & Ham 



mer 



(2010) to investigate whether the Large 



Magellanic Cloud could be related to some 
ejected material coming from M31. 

Note added into proofs: The specific frequency 
of globular clusters in M31 appears to be three or 
four times greater than it is in the Galaxy [van den| 
Bergh ( 2010 ), i.e., it is expected if a major merger 



occurred in the first galaxy. The fact that the glob- 
ulars in M31 and in the Galaxy have almost ex- 
actly the same half light radii supports the notion 
that M31 formed by merger of two Galaxy-sized 
ancestors (Sidney van den Bergh, private commu- 
nication) . 
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